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' We propose a new set of measurements which can be performed at the LHC using roman pot 

detectors. The method exploits excitation curves in central diffractive pair production, and is illus- 

. trated using the examples of the W boson and top quark mass measurements. Further applications 

, are mentioned. 
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I. INTRODUCTION 



We propose a new method to measure heavy particle properties via double photon and double pomeron exchange 
(DPE), at the LHC. In this category of events, the heavy objects are produced in pairs, whereas the beam particles 
often leave the interaction region intact, and can be measured using very forward detectors. 

If the events are exclusive, i.e., if no other particles are produced in addition to the pair of heavy objects and the 
1 outgoing protons, the proton measurement gives access to the photon-photon or pomeron-pomeron centre-of-mass [1], 
and the dynamics of the hard process can be accurately studied. In particular, one can observe the threshold excitation 

■ and attempt to extract the heavy particle's mass, or study the particle's (possibly energy-dependent) couplings by 
measuring cross-sections and angular distributions. As examples of this approach, we give a detailed account of the 

• W boson and top quark mass measurement at production threshold. The method can easily be extended to other 
Q_i' heavy objects. 

We also consider inclusive double pomeron exchange, i.e., events where other particles accompany the heavy system. 
This production mode is to be expected in central diffractive tt production, since inclusive double diffractive dijet 
production has effectively been observed at the Tevatron [2] . The cross-section measurement of [2] allows for a rough 

■ estimate of the LHC expectation [3] . The interest of such events is reviewed. 
The paper is organised as follows. We start by giving the theoretical formulation of W + W~ production (via QED) 

and of tt production (in both exclusive and inclusive DPE). We then describe the event generation, the simulation of 
detector effects, and the cuts used in the analysis. The following part of the paper describes in detail the threshold 
scan method, in a twofold version ( "turn-on" and "histogram" fits), and its application to the W boson and top quark 
mass measurements. We then conclude on the method in general, on the above mass measurements in particular, and 
mention a number of further applications. 



II. THEORETICAL FORMULATION 



Pair production of W bosons and top quarks in QED and double pomeron exchange are described in detail in 
this section. W + W~ pairs are produced in photon-mediated processes, which are exactly calculable in QED. There 
is basically no uncertainty concerning the possibility of measuring these processes at the LHC. On the contrary, tt 
events, produced in exclusive double pomeron exchange, suffer from theoretical uncertainties since exclusive diffractive 
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production is still to be observed at the Tevatron, and other models lead to different cross sections, and thus to a 
different potential for the top quark mass measurement. However, since the exclusive kinematics are simple, the model 
dependence will be essentially reflected by a factor in the effective luminosity for such events. 

By contrast, the existence of inclusive double pomeron exchange — in other words, when the pomeron remnants 
carry a part of the available center-of-mass energy — is certain since it has been observed already in experiments. We 
will mention at the end how these events could be used and the interest of their experimental determination. We will 
briefly analyse their impact on the tt threshold scan but we postpone a precise study of such events to a forthcoming 
publication. 
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FIG. 1: W + W~ (QED) and tt (QCD) exclusive production. Left: double photon exchange process. Right: exclusive double 
pomeron exchange in the Bialas-Landshoff model via "gluons in the pomeron"; The grey band represents the rapidity gap 
survival suppression factor (see text). 



A. W + W production via double photon exchange 

The QED process rates are obtained from the following cross section formula 

^ CT (pp^ p w+w-p) = ^77 -> w+w- dnj dnl , 



where the Born 77 — > W + W cross-section reads [4] 
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where Mww is the total W + W mass. The photon fluxes dn 1 are given by [5] 
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in the usual dipole approximation for the proton electromagnetic form factors. u> is the photon energy in the laboratory 
frame, q 2 the modulus of its mass squared in the range 



r 2 2 



- 2 q 2 1 = 

dmaxj — 



E(E-cu) ' q 2 



(6) 



where E and m are the energy and mass of the incident particle and t max = (m^ v /M^ w ) maa . is defined by the 
experimental conditions. 

The QED cross section da(j>p — > p W + W~p) is a theoretically clear prediction. One should take into account 
however, two sources of (probably mild) correction factors. One is due to the soft QCD initial state radiation between 
incident protons which could destroy the large rapidity gap of the QED process. It is present but much less pronounced 
than for the rapidity gap survival for a QCD hard process (see the discussion in the next subsections), thanks to the 
large impact parameter implied by the QED scattering. The second factor is the QCD gg — > W + W~ exclusive 
production via higher order diagrams. This has been evaluated recently [6] for standard (non diffractive) production 
to give a 5% correction factor. The similar calculation for the diffractive W + W~ production by comparison with the 
QED process is outside the scope of our paper but deserves to be studied together with the "inclusive" background 
(W + W~+hadrons) it could generate. 



B. Exclusive diffractive production of tt events 

Let us introduce the model [7, 8] we shall use for describing exclusive tt production in double diffractive production. 
This process is depicted in Fig. I (right). 

In [7], the diffractive mechanism is based on two-gluon exchange between the two incoming protons. The soft 
pomeron is seen as a pair of gluons non-perturbatively coupled to the proton. One of the gluons is then coupled 
perturbatively to the hard process (either the Higgs boson, or the tt pair, see Fig. 1), while the other one plays the 
role of a soft screening of colour, allowing for diffraction to occur. The corresponding cross-sections for qq and Higgs 
boson production read: 




datr(s) = qr ( ) <5 (2) [ 5>« + hi) | J] d 2 v t d 2 hdm £ Q '" ? exp(-2Att«?) *tt , (7) 

»=1,2 

where the variables Vi and h respectively denote the transverse momenta of the outgoing protons and of top quarks, 
£i are the proton fractional momentum losses, and rn are the quark rapidities, 

_ tt da _ p(l - p) _ Am 2 
u ~ 24 dt ~ m 2 T1 m 2 T2 ' 9 ~ [ ' 

is the hard gg — > tt cross-section. 

In the model, the soft pomeron trajectory is taken from the standard Donnachie-Landshoff parametrisation [9], 
namely a(t) = 1 + e + a't, with e w 0.08 and a' w 0.25GcV~ 2 . X t i and the normalization C^ c are kept as in the 
original paper [7]. Note that, in this model, the strong (non perturbative) coupling constant is fixed to a reference 
value G 2 /An — I, reflecting the lack of knowledge of the absolute normalization of exclusive DPE processes. 



C. Inclusive diffractive production of tt events 

It is convenient to introduce also the model for central inclusive diffractive production [3] applied to tt dijets. One 
writes 

dair 1 = C tl (^t) I £ Vi + kA Hld^dxfdVh^h^i^X^)} a tl G P (x{ , y)G P {^ 2 , y) . (9) 

^ tt / \i=l,2 J i=l,2 

In the above, x\ are the pomeron's momentum fractions carried by the gluons or quarks involved in the hard pro- 
cess, and the Gp is the gluon energy density in the pomeron, i.e., the gluon density multiplied by xf. We use as 
parameterisations of the pomeron structure functions the fits to the diffractive HERA data performed in [10]. The 
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normalization C t i is obtained from the description [3] of the jet-jet diffractive cross-section at the Tevatron [2]. The 
hard cross-section a t i to be considered is now 




to be distinguished from a t i (8) due to inclusive characteristics [3] . 

D. Rapidity Gap Survival 

In order to select exclusive diffractive states, such as for W + W~ (QED) and tt (exclusive, QCD), it is required to 
take into account the corrections from soft hadronic scattering. Indeed, the soft scattering between incident particles 
tends to mask the genuine hard diffractive interactions at hadronic colliders. The formulation of this correction [11] 
to the scattering amplitude -4(ww,tt) consists in considering a gap survival probability (SP) function S such that 

A(p T i,PT2, A$) = {1 + Asp} x.4(ww,tt) =5x^ (WWitt -) = J d 2 k T «S(k T ) -4( W w,tt)(PTi-k T , p T 2 + k T ) , (11) 

where pn,2 are the transverse momenta of the outgoing p,p and A<I> their azimuthal angle separation. Asp is the 
soft scattering amplitude. 

The correction for the QED process is present but much less pronounced than for the rapidity gap survival for a 
QCD hard process, thanks to the large impact parameter implied by the QED scattering. In a specific model [12] the 
correction factor has been evaluated to be of order 0.9 at the LHC for 77 — > H. It is evaluated to be of order 0.03 for 
the QCD exclusive diffractive processes at the LHC. 

III. EXPERIMENTAL CONTEXT 

A. The DPEMC Monte Carlo 

A recently developed Monte-Carlo program, DPEMC [13], provides an implementation of the W + W~ and tt events 
described above in the QED and both exclusive and inclusive double pomeron exchange modes. It uses HERWIG [14] as 
a cross-section library of hard QCD processes, and when required, convolutes them with the relevant pomeron fluxes 
and parton densities. The survival probabilities discussed in the previous section (respectively 0.9 for double photon 
and 0.03 for double pomeron exchange processes) have been introduced at the generator level. The cross section at 
the generator level for W + W~ QED and exclusive diffractive tt production is found to be 55.9 fb and 40.1 fb for a 
mw mass of 80.42 GeV and a top mass of 174.3 GeV after applying the survival probabilities. 

B. Roman pot detector positions and resolutions 

A possible experimental setup for forward proton detection is described in detail in [15]. We will only describe its 
main features here and discuss its relevance for the W boson and top quark masses measurements. 

In exclusive DPE or QED processes, the mass of the central heavy object can be reconstructed using the roman 
pot detectors and tagging both protons in the final state at the LHC. It is given by M 2 — Ci&s, where & are the 
proton fractional momentum losses, and s the total center-of-mass energy squared. In order to reconstruct objects 
with masses in the 160 GeV range (for W + W~ events) in this way, the acceptance should be large down to £ values 
as low as a few 10 -3 . For tt events, an acceptance down to 10 -2 is needed. The missing mass resolution directly 
depends on the resolution on £, and should not exceed a few percent to obtain a good mass resolution. 

These goals can be achieved if one assumes two detector stations, located at ~ 210 m, and ~ 420 m [15] from the 
interaction point 1 . The £ acceptance and resolution have been derived for each device using a complete simulation of 
the LHC beam parameters. The combined £ acceptance is close to ~ 60% at low masses (at about twice mw), and 



A third position at 308 m is often considered as well but is more difficult from a technological point of view at the LHC and was not 
considered for this study. 
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90% at higher masses starting at about 220 GeV. for £ ranging from 0.002 to 0.1. The acceptance limit of the device 
closest to the interaction point is £ > £ TOi „ =0.02. Let us note also that the acceptance for tt events goes down to 
20% if only roman pots at 210 m are present since most of the events are asymmetric (one tag at 420m and another 
one at 210m). 

Our analysis does not assume any particular value for the £ resolution. We will discuss in the following how the 
resolution on the W boson or the top quark masses depend on the detector resolutions, or in other words, the missing 
mass resolution. 



C. Experimental cuts 

This section summarises the cuts applied in the remaining part of the analysis. As said before, both diffracted 
protons are required to be detected in roman pot detectors. 

The triggers which will be used for the W + W~ and tt events will be the usual ones at the LHC requiring in addition 
a positive tagging in the roman pot detectors. 

The experimental offline cuts and their efficiencies have been obtained using a fast simulation of the CMS detector 
[16] as an example, the fast simulation of the ATLAS detector leading to the same results. If we require at least one 
lepton (electron or muon) with a transverse momentum greater than 20 GeV and one (two) jet with a transverse 
momentum greater than 20 GeV (40 GeV) for W + W~ (tt) to be reconstructed in the acceptance of the main detector 
in addition to the tagged protons 2 , we get an efficiency of about 50% for tt events, and 30% for W + W~ events. We 
give the mass resolution as a function of luminosity in the following after taking into account these efficiencies. If the 
efficiencies are found to be higher, the luminosities have to be rescaled by this amount. 

IV. THRESHOLD SCAN METHODS 

A. Explanation of the methods 

We study two different methods to reconstruct the mass of heavy objects double diffractively produced at the 
LHC. As we mentioned before, the method is based on a fit to the turn-on point of the missing mass distribution at 
threshold. 

One proposed method (the "histogram" method) corresponds to the comparison of the mass distribution in data 
with some reference distributions following a Monte Carlo simulation of the detector with different input masses 
corresponding to the data luminosity. As an example, we can produce a data sample for 100 fb _1 with a top mass of 
174 GeV, and a few MC samples corresponding to top masses between 150 and 200 GeV by steps of. For each Monte 
Carlo sample, a x 2 value corresponding to the population difference in each bin between data and MC is computed. 
The mass point where the % 2 is minimum corresponds to the mass of the produced object in data. This method has 
the advantage of being easy but requires a good simulation of the detector. 

The other proposed method (the "turn-on fit" method) is less sensitive to the MC simulation of the detectors. 
As mentioned earlier, the threshold scan is directly sensitive to the mass of the diffractively produced object (in the 
W + W~ case for instance, it is sensitive to twice the W mass). The idea is thus to fit the turn-on point of the missing 
mass distribution which leads directly to the mass of the produced object, the W boson. Due to its robustness, this 
method is considered as the "default" one in the following. 

To illustrate the principle of these methods and their achievements, we apply them to the W boson and the top 
quark mass measurements in the following, and present in detail the reaches at the LHC. They can be applied to 
other threshold scans as well. 



B. W mass measurement 

In this section, we will first describe the result of the "turn-on fit" method to measure the W mass. As we mentioned 
earlier, the advantage of the W + W~ processes is that they do not suffer from any theoretical uncertainties since this 



2 The double pomeron exchange background to the signal is found to be small, and will more correspond to misidentifications of jets as 
leptons in the main detector. Since this is difficult to evaluate precisely using a fast simulation of the detector, and this is quite small 
compared to the signal, we decided not to incorporate it in the following study. 
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FIG. 2: Two examples of fits to missing mass reference distributions with a resolution of the roman pot detectors of 1 GeV 
(left) and 3 GeV (right). We see on these plots the principle and the accuracy of the "turn-on fits" to the MC at threshold. 




FIG. 3: Calibration curves (see text) for two different roman pot resolutions of 1 GeV (left) and 3 GeV (right). We notice that 
the calibration can be fitted to a linear function with good accuracy. The dashed line indicates the first diagonal to show the 
shift clearly. 



is a QED process. The W mass can be extracted by fitting a 4-parameter ! turn-on' curve to the threshold of the mass 
distribution (c.f. Ref. [17]): 



T = P X 



— +1 



P 4 • (12) 



Pi is the amplitude, P2 the inflexion point, P3 the width of the turn-on curve, and P4 is a vertical offset, x being 
the missing mass. With a detector of perfect resolution, P2 would be equal to twice the W mass. However, the finite 
roman pot resolution leads to a shift between P2 and 2mw which has to be established using a MC simulation of the 
detector for different values of its resolution. This shift is only related to the method itself and does not correspond to 
any error in data. For each value of the W input mass in MC, one has to obtain the shift between the reconstructed 
mass (P2/2) and the input mass, which we call in the following the calibration curve. It is assumed for simplicity that 
P2 is a linear function of mwi which is a good approximation as we will see next. In order to determine the linear 
dependence between P2 and mw, calibration curves are calculated for several assumed resolutions of the roman pot 
detectors. The calibration points are obtained by fitting T to the mass distribution of high statistics samples (100 000 
events) for several values of mw. An example is given in Fig. 2 for two resolutions of the roman pot detectors. The 
difference between the fitted values of P2/2 and the input W masses are plotted as a function of the input W mass and 
are then fitted with a linear function. To minimise the errors on the slope and offset, the difference P2/2 — 80.42 GeV 
is plotted versus mw (Fig- 3). 

To evaluate the statistical uncertainty due to the method itself, we perform the fits with some 100 different "data" 
ensembles. For each ensemble, one obtains a different reconstructed W mass, the dispersion corresponding only to 
statistical effects. The expected statistical uncertainty on the actual measurement of the W mass in data is thus 
estimated with these ensemble tests for several integrated luminosities and roman pot resolutions. Each ensemble 
contains a number of events that corresponds to the expected event yield for a given integrated luminosity, taking into 
account selection and acceptance efficiencies. The turn-on function T is fitted to each ensemble. Only the parameters 
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FIG. 4: Distribution of the fitted value of the W mass from ensemble tests. Left: corresponding to 150 fb 1 , right: corresponding 
to 300 fb~ . We note the resolution obtained on the W mass for these two luminosities. 
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FIG. 5: Expected statistical uncertainty on the W mass as a function of luminosity for three different roman pot resolutions. 

Pi and P2 are allowed to float, P3 and Pa are fixed to the average values obtained from the fits for the calibration 
points. 

In order to obtain the fitted estimate for the W mass, m^, in each ensemble, the fit value of P2 is corrected with the 
calibration curve that corresponds to the roman pot resolution. For each resolution is histogrammed as shown 
in Fig. 4. The distributions are fitted with a Gauss function where the width corresponds to the expected statistical 
uncertainty of the W mass measurement. Fig. 5 shows the expected precision as a function of the integrated luminosity 
for several roman pot resolutions. With 150 fb -1 the expected statistical uncertainty on mw is about 0.65 GeV when 
a resolution of the roman pot detectors of 1 GeV can be reached. With 300 fb -1 the expected uncertainty on mw 
decreases to about 0.3 GeV. 

We notice of course that this method is not competitive to get a precise measurement of the W mass, which would 
require a resolution to be better than 30 MeV. However, this method can be used to align precisely the roman pot 
detectors for further measurements. A precision of 1 GeV (0.3 GeV) on the W mass leads directly to a relative 
resolution of 1.2% (0.4%) on £ using the missing mass method. This calibration will be needed, for instance, to 
measure the top mass as proposed in the next section. 

Let us now present the result on the "histogram" method, which is an alternative approach to determine the W 
mass. The same high statistics templates used to derive the calibration curves are fitted directly to each ensemble 
(see Fig. 6 left). The \ 2 is defined using the approximation of poissonian errors as given in Ref. [18]. Each ensemble 
thus gives a x 2 curve which in the region of the minimum is fitted with a fourth-order polynomial (Fig. 6 right). The 
position of the minimum of the polynomial, m w m , gives the best value of the W mass and the uncertainty cr(mw) is 
obtained from the values where \ 2 = Xmin + 1- The mean value of cr(mw) for ah ensembles are quoted as expected 
statistical uncertainties. 

The expected statistical errors on the W mass using histogram fitting are comparable to those using the function 
fitting method. However, the former turns out to be more sensitive to the resolution of the roman pot detectors. 

C. Top mass measurement 

The method to extract the top mass is the same as for the W mass described in the previous section. The theoretical 
cross section is not as well known as for the W and is model dependent. Our study assumes the Bialas Landshoff 
model for exclusive tt production. For tt events the width of the turn-on curve is considerably larger than for WW 
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total pot mass(GeV) m w (GeV) 

FIG. 6: Left: Example of the histogram-fitting method. We see the difference between the "data" sample (full histogram 
with error bars, mw = 80.42 GeV) and a reference histogram (dashed line, mw = 85.42 GeV). Right: Example of the \ 2 
distribution in one ensemble. 
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FIG. 7: Expected statistical precision of the W mass as a function of the integrated luminosity for various resolutions of the 
roman pot detectors using the histogram-fitting method. 

events (Fig. 8, left), resulting in a larger offset between the actual turn on and the inflexion point of the fit function 3 . 
The calibration curve for a resolution of the roman pot detectors of 1 GeV is displayed in Fig. 8, right. 

Ensemble tests for integrated luminosities of 50, 75, 100 and 200 fb" 1 and roman pot detector resolutions of 1 GeV, 

2 GeV and 3 GeV yield the results shown in Fig. 9, left. Resolutions of the roman pot detectors between 1 GeV and 

3 GeV give similar statistical uncertainties on the top quark mass which is due to the fact that the main limiting 
effect on resolution is statistics. With 100 fb" 1 the expected statistical precision is about 1.6 GeV and gets improved 
to about 0.65 GeV with 300 fb" 1 . 

The results have also been cross-checked using the histogram fitting method which was found to yield very similar 
expected uncertainties as the function fitting method (Fig. 9, right). 

V. OUTLOOK AND PROSPECTS 

In this section, we discuss other applications of the threshold scan method. Detailed analysis is postponed to 
forthcoming papers [19]. 

As we mentioned before, the cross section of exclusive top pair production at the LHC is still uncertain, and 
predictions will be constrained by the incoming results form the Tevatron, especially from the D0 experiment where 
it is possible to detect double tagged events. On the contrary, inclusive double pomeron exchange has already been 
observed, and top quark pair production in this mode is fairly certain at the LHC. In this case, the threshold excitation 
is sensitive to quark and gluon densities at high pomeron momentum fraction, so that these events provide a rather 
unique opportunity to study structure functions near the endpoint. 



3 Note in addition that the top quark width is not included in Herwig and thus in our study. However, this effect is expected to be small. 
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FIG. 8: Fit to a reference mass distribution with mt = 175 GeV (left) and calibration curve for a roman pot resolution of 
1 GeV, the diagonal is displayed in dashed line to show the difference (right). Note the larger difference between the calibration 
curve and the diagonal compared to the case of W + W~ production. 




400 
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FIG. 9: Expected statistical precision of the top mass as a function of the integrated luminosity for various resolutions of the 
roman pot detectors (full line: resolution of 1 GeV, dashed line: 2 GeV, dotted line: 3 GeV). Left: function fitting method, 
right: histogram fitting method (the three curves for different roman pot resolutions lead to the same results and cannot be 
distinguished in the figure). 

To illustrate this point, we give in Fig. 10 the missing mass distributions at the generator level using the DPEMC 
Monte Carlo for the exclusive tt events (full line) and the results on the Bialas-Landshoff inclusive tt production for 
two different gluon densities in the pomeron (dashed line: fit 1, dotted line: fit 2, see Ref. [20]). Fit 2 in Ref [20] 
leads to a more prominent gluon at high (3 than fit 1. We see that the missing mass distribution is directly sensitive 
to the parton distributions in the pomeron. In Fig. 11, we display the differences between the exclusive tt events in 
full line and the result of the factorisable POMWIG model (dotted line), and the non factorisable one based on the 
Bialas-Landshoff approach. We see again that the missing mass distribution, and thus the threshold analysis can help 
distinguishing between the models. 

Another application of exclusive pair-production consists in measuring the mass of stops and sbottoms, provided 
these particles exist and can be produced in pairs at the LHC. 

Finally, W pair-production in central diffraction gives access to the couplings of gauge bosons. Namely, as we 
mentioned already, W + W~ production in two-photon exchange is robustly predicted within the Standard Model. 
Any anomalous coupling between the photon and the W will reveal itself in a modification of the production cross 
section, or by different angular distributions. Since the cross-section of this process is proportional to the fourth power 
of photon-W coupling, good sensitivity is expected. 

VI. CONCLUSION 

Recent work on DPE has essentially focused on the Higgs boson search in the exclusive channel. In view of the 
difficulties and uncertainties affecting this search [8] , we highlight new aspects of double diffraction which complement 
the diffractive program at the LHC. 
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FIG. 10: Missing mass distributions at the generator level using the DPEMC Monte Carlo for the exclusive tt events in full 
line and the results on inclusive tt production for two different gluon densities in the pomeron (dashed line: fit 1, dotted line: 
fit 2, see Ref. [20]). 



3500 
3000 
2500 
2000 
1500 
1000 
500 



400 600 800 1000 1200 1400 

Central mass (GeV) 



FIG. 11: Missing mass distributions at the generator level using the DPEMC Monte Carlo for the exclusive tt events in full 
line, the inclusive events from model [3, 8] (dashed line), and [21] (dotted line), see text. 

In particular, QED W pair production provides a certain source of interesting diffractive events. Inclusive tt 
production via double pomeron exchange is also an open channel and will provide interesting information on a poorly 
known aspect of structure functions. These robust channels will help and accompany the understanding of the more 
intriguing and challenging problem of exclusive double diffraction. 

In this paper, we have advocated the interest of threshold scans in double pomeron exchange. This method 
considerably extends the physics program at the LHC. To illustrate its possibilities, we described in detail the W 
boson and the top quark mass measurements. The precision of the W mass measurement is not competitive with 
other methods, but provides a very precise calibration of the roman pot detectors. The precision of the top mass 
measurement is however competitive, with an expected precision better than 1 GeV at high luminosity. Other 
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promising applications remain to be investigated. 
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